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ABSTRACT 

Magnetic fields in upper main-sequence stars, white dwarfs, and neutron stars are 
known to persist for timescales comparable to their lifetimes. From a theoretical per- 
spective this is problematic, as it can be shown that simple magnetic field configu- 
rations are always unstable. In non-barotropic stars, stable stratification allows for a 
much wider range of magnetic field structures than in barotropic stars, and helps sta- 
bilize them by making it harder to induce radial displacements. Recent simulations by 
Braithwaite and collaborators have shown that, in stably stratified stars, random ini- 
tial magnetic fields evolve into nearly axisymmetric configurations with both poloidal 
and toroidal components, which then remain stable for some time. It is desirable to 
provide an analytic study of the stability of such fields. We write an explicit expres- 
sion for a plausible equilibrium structure of an axially symmetric magnetic field with 
both poloidal and toroidal components of adjustable strengths, in a non-barotropic 
static fluid star, and study its stability using the energy principle. We construct a 
displacement field that should be a reasonable approximation to the most unstable 
mode of a toroidal field, and confirm Braithwaite's result that a given toroidal field 
can be stabilized by a poloidal field containing much less energy than the former, as 
given through the condition -Epoi/£^tor ^ 2a£'tor/£'grav, where E-po\ and i?tor are the 
energies of the poloidal and toroidal fields, respectively, and -Egrav is the gravitational 
binding energy of the star. We find that a « 7.4 for main-sequence stars, and a ~ 200 
for neutron stars. Since i^poiZ-E-grav ^ 1, we conclude that the energy of the toroidal 
field can be substantially larger than that of the poloidal field, which is consistent 
with the speculation that the toroidal field is the main reservoir powering magnetar 
activity. 

Key words: instabilities - magnetic fields - MHD - stars: magnetic field - stars: 
neutron - white dwarfs. 



1 INTRODUCTION 

Upper main-sequence stars, white dwarfs, and neutron stars are known to possess magnetic fields that persist for long 
periods of time, comparable to their lifetimes. Since convection does not play an important role in these objects, dynamo 
generation of magnetic fields is not expected during most of their lives. As a consequence, their magnetic fields must be 
in stable hydromagnetic equilibrium. However, from a theoretical perspective this poses a problem, since it can be shown 
that simple magnetic field configurations consisting of purely poloidal (meridional) or purely toroidal (azimuthal) fields are 
always unstable. In particular, Tayler (1973) showed that toroidal fields are prone to the interchange (axisymmetric) and kink 
(non-axisymmetric) instabilities. Markey & Tayler (1973) and Wright (1973) showed similarly that purely poloidal fields, with 
some field lines closing inside the star, are also unstable near the neutral line, where the poloidal field vanishes. Flowers & 
Ruderman (1977) discussed another large-scale instability of poloidal fields, illustrated by the fact that, when two magnets 
are aligned, they will tend to orient in opposite direction to one another. Therefore, a rotation of an entire hemisphere of 
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a star, cut along a plane containing the axis, should lead to the lowering of the overall energy, as we have demonstrated 
mathematically (Marchant, Reisenegger & Akgiin 2011). 

On the other hand, even in the most strongly magnetized stars, the magnetic (Lorentz) force inferred from the surface 
field strengths is still typically a million times weaker than the hydrostatic force due to pressure and gravity. Therefore, a 
small perturbation in the non-magnetic background equilibrium could be sufficient to balance the magnetic force. In the 
radiative envelopes of massive stars and in the interiors of degenerate stars, matter is non-barotropic (i.e. pressure depends 
on a second quantity, such as chemical composition or specific entropy, in addition to density) and stably stratified, allowing 
for a wider range of magnetic field structures than found in barotropic fluids (i.e. those where pressure can be expressed as a 
function of density only) (Reisenegger 2009). Stable stratification also helps stabilize the magnetic field by making it harder 
to induce radial displacements of the fluid. This effect was included by Tayler (1973); however, by itself it is not sufficient to 
completely stabilize a purely toroidal (or purely poloidal) magnetic field. Recent simulations for stably stratified stars have 
demonstrated that initially random magnetic fields tend to evolve into nearly axisymmetric configurations with both poloidal 
and toroidal components of comparable strength, which then remain stable for several Alfven times (Braithwaite & Spruit 
2004; Braithwaite & Nordlund 2006; Braithwaite 2009). Our goal is to provide an analytic justification for the stability of 
such fields, and to understand how the poloidal and toroidal components can help stabilize each other. 

Before considering the stability, we must first determine the equilibrium structure of the magnetic field. In barotropic 
stars, the equilibrium form of the magnetic field is severely restricted and is given as the solution of a differential equation (the 
so-called Grad-Shafranov equation, as discussed, for example, in Chandrasekhar & Fermi 1953; Ferraro 1954; Liist & Schliiter 
1954; Prendergast 1956; a detailed discussion is also given in Akgiin & Wasserman 2008). On the other hand, in realistic, 
non-barotropic stars, the hydrostatic force includes a buoyancy term, which acts as a restoring force for stably stratified 
fluids. Thus, the only restriction that remains for an axisymmetric field in a non-barotropic fluid is that the magnetic force 
cannot have an azimuthal (cf)) component, since no counterpart exists in the hydrostatic force that can act to balance it. In 
addition, the equilibrium magnetic field needs to satisfy boundary conditions at the surface and regularity conditions at the 
center of the star. We can construct simple polynomial forms for the scalar functions that describe the poloidal and toroidal 
components of the magnetic field, consistent with these requirements. 

Once we know the equilibrium structure of the magnetic field, we can examine its stability, for which we use the energy 
principle developed by Bernstein et al. (1958). In this method, one considers the energy of perturbations around the magnetic 
equilibrium. If this energy is positive, then the equilibrium is stable, and vice versa. We then consider the problem of 
constructing a displacement field that gives rise to instabilities in a purely toroidal field configuration. The hydrostatic and 
toroidal parts of the energy can be examined analytically for stability, and can be minimized with respect to the azimuthal 
component of the displacement field, in an analogous manner to Tayler (1973). Once we have found this minimum, we add 
the poloidal part of the energy and determine how strong the poloidal field must be in comparison to the toroidal field in 
order to stabilize the field. 

The outline of this paper is as follows. In iJ2]we discuss the equilibrium structure of the star and the magnetic field. We 
first construct sample equilibrium profiles for the pressure, density, and gravitational potential, which, while being sufficiently 
simple, have all the desirable qualities. Next, we consider the structure of the poloidal and toroidal fields and discuss their 
properties. We then construct a simple magnetic field that satisfies the boundary and regularity conditions. In [J2]we consider 
the stability of the magnetic field thus constructed using the energy principle approach. We calculate the contributions to 
the energy due to the fluid, and due to the poloidal and toroidal components of the magnetic fleld. We give a proof that all 
physically relevant, purely toroidal fields are unstable. We discuss the implications of stable stratification on the displacement 
fleld. We construct a particular displacement fleld that makes the sum of the hydrostatic and toroidal parts of the energy 
negative, yielding an instability, and then show how the addition of a poloidal field eliminates this instability. In !2]we present 
our conclusions. 



2 EQUILIBRIUM 

In realistic stars, the stress due to the magnetic field is much weaker than the hydrostatic terms due to pressure and gravity 
(e.g. Reisenegger 2009). The background equilibrium in the absence of magnetic fields and rotation is spherically symmetric 
and is given by Euler's equation, 

VPo + f?oV$o = , (1) 

where P is pressure, g is density, and <I> is gravitational potential. Throughout this paper, we will denote the spherically 
symmetric non-magnetic background quantities with the subscript 0. The gravitational potential is given in terms of the 
density by Poisson's equation, 

V^c&o = 4TvGgo . (2) 
The magnetic fleld B changes the background quantities slightly, and the new equilibrium is given by 

VP + gV$ = ^^-^ , (3) 

c 

where J — cV x B/Att is the current density. We can express the small changes due to the magnetic fleld as Eulerian 
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perturbations, and write P = Po + Pi, and similarly for g and "1>. Then, we can rewrite the above equation, working to first 
order in the perturbations, as 

J X B 

VPl + £>lV$o + £>oV$i = — - — . (4) 

In the often used, idealized assumption of barotropic fluids, there is a unique relation between pressure and density, 
which holds throughout the application of small perturbations. Therefore, we can write the pressure as a function of the 
density. This allows us to express the left-hand side of equation ([Sjl (and consequently equation [Jjl as a gradient of the form 
VP-I- (jV<1? — pV(-ff + $), where dH{g) — dP^g) / g. This implies that V x ( J x B/ gc) = 0, so the magnetic acceleration must 
also be expressible as a gradient. This is a strong constraint and greatly restricts the possible choice of the magnetic field 
in equilibrium (Chandrasekhar & Fermi 1953; Ferraro 1954; Liist & Schliiter 1954; Prendergast 1956; Akgiin & Wasserman 
2008; HaskeU et al. 2008). 

On the other hand, in non-barotropic fluids, pressure depends on at least one additional quantity, as well as density. In 
white dwarfs and in the radiative zones of non-degenerate stars, the dominant additional quantity is the speciflc entropy, and 
in neutron stars it is the composition (fraction of protons or other "impurities"; Reisenegger 2009). For long equilibration 
times, any changes induced in the background quantities will imply that a simple relation between pressure and density no 
longer exists. Consequently, the left-hand sides of equations Q and Q are not expressible as gradients. Therefore, unlike the 
barotropic case, we do not require that the magnetic acceleration be expressible as a gradient. Instead, the only constraint for 
axisymmetric fields is the much less restrictive requirement that the 4> component of the magnetic force density vanish, since 
there is no such component in the hydrostatic part that can balance it (Chandrasekhar & Prendergast 1956; Mestel 1956). 



2.1 Non-magnetic equilibrium 

In this section, we will consider a simple model for the non-magnetic background equilibrium quantities. The derivations that 
follow in the subsequent sections do not rely on the speciflc model, but it will be needed later in the calculation of numerical 
estimates. A density profile that is simple enough and does not deviate by more than a few percent from an ti = 1 polytrope 
(Mastrano et al. 2011), is 

go{x) = gc{l - x^) , (5) 

where g^ is the central density, and we define a dimensionless radial coordinate x = r/Ri,, where P* is the stellar radius. The 
mass enclosed within radius x is given by 

moix) = 47rP^ T go{x)x''dx = l^!-?*^(5a;^ - 2.x^) . (6) 
Jo 15 

If the total mass of the star is denoted by Af*, then the central density is gc = 15M*/87rPj. The gravitational potential inside 
the star is given by Poisson's equation (equation [21) , 



Here, the gravitational potential at the center is chosen to be zero. From Euler's equation (equation [T]), we find that the 
pressure is given by 

P„(,) = P_£^^M^:g^rf,.P,(^i_^ + 2x*-^) . (8) 

The value of the central pressure Pc is determined by requiring that at the surface Po (1) ~ 0, which yields P^ — AnGg^R't /15 = 
15GM^/167rP*. The profiles of the background quantities Po, go and $0 are shown in Fig.[T] The gravitational binding energy 
of the star is 

-Bgrav = 47rGP* / xgo(x)mo[x)dx = —— — = — . (9) 

.In iK* ^1 



2.2 Magnetic field structure 

The magnetic field is divergenceless, therefore quite generally it can be expressed as the sum of a poloidal and a toroidal 
component, each completely described by a single scalar function (Chandrasekhar 1981). These functions are analogous to 
the stream functions describing incompressible flows in hydrodynamics. An axisymmetric magnetic field can be written in 
spherical coordinates (r, 9, (j>) as 

B = Bpoi + Btor = Va(r, 61) x V(?i + /3(r, 6») V0 . (10) 
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Here, we make use of the relation — (j) / {r sin 6) , which simplifies the calculation of curls. The current densities corre- 
sponding to each component are given by 

= V X Bpo, = -AaV<^ , 
= V X Btor = V/3 X V<^ . 

c 

We have introduced the so-called Grad-Shafranov operator, defining the cylindrical radius as = rsinO, 

A = r^'V-{z.-'V)^d'^ + '^do(^) . (12) 

The curl of a poloidal field is a toroidal field, and the curl of a toroidal field is a poloidal field. Therefore, Jpoi is actually a 
toroidal field, and Jtor is a poloidal field. 

We have Jpoi || -Btor || therefore the term J-po\ x -Btor always vanishes in the Lorentz force. On the other hand, for 
the poloidal components we have Jtor -L (f> and Bpoi -L 4>, which implies that Jpoi x Bpoi -L 4> and Jtor x i?tor -L 4>, while 
Jtor X Bpoi II 4>- However, in axisymmetric equilibrium the Lorentz force cannot have a (f> component, as implied by equation 
((3}. Therefore, we must also have Jtor || Spoi, or equivalently Vq || V/3, which implies that /3 can be expressed as a function 
of a, i.e. P = /3(a). Thus, the Lorentz force can be written as the sum of a term entirely due to the poloidal field, and one 
entirely due to the toroidal field, /„j^g — f^^i + /ton where 



47i-/poi = (V X Bpoi) X Bpoi = -vj ^AaVa 



47r/tor = (V X Btor) X Btor = -Vj''^ p = -tlJ-^/S^Va 

da 

Note that the two terms are poloidal and parallel. Moreover, they are perpendicular to the magnetic surfaces (defined as the 
surfaces of constant a and /3, which contain all the field lines). 



(13) 



2.3 Poloidal field 

In this section, we derive a simple profile for an axisymmetric poloidal magnetic field that conforms to certain boundary and 
regularity conditions. In particular, we impose that there are no surface currents (which would be dissipated very quickly), 
implying that the poloidal field is continuous across the surface. We assume that the current density drops continuously 
towards the surface, as the number density of charged particles should be decreasing with the mass density. Moreover, the 
magnetic field and current density should remain finite and continuous everywhere in the interior, and in particular at the 
center of the star. We can then construct a polynomial solution for the scalar function that describes the poloidal magnetic 
field, consistent with these requirements. 

Writing the dimensional part of the magnetic field explicitly in terms of some constant Bo, the poloidal field can be 
expressed as Bpoi = BoVa x (equation llOf) . Here, hats denote that the operators are with respect to the dimensionless 
radial coordinate x = r/7?*, and & is also dimensionless. We assume that the field outside the star is that of a point dipole. 
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which is of the form Bdip oc x ^(2rcos6' + Osind), corresponding to a{x,6) oc sin^6'/a;. In order to match the angular 
dependence of the field on the surface, we choose 

a{x,e) = f{x)sm''e . (14) 



The poloidal field becomes (equation [10} 

-Bpoi = -BoVa X = Bo 
The current density is (equation I lip 



2 fix) cose ^ _ /(x)sin6' ^' 



(15) 



where, from equation (|12|) . we have 



^^ = VxBpo, = -|^AqV,/>, (16) 
c K* 



Aa = ( /" - |{ ) sin^e . (17) 



Outside the surface, the current density is zero, which implies that 

/" = 1^ for x>l. (18) 

Plugging in a trial solution of the form f (x x" , we find that the solutions are s = — 1 and s — 2. Thus, outside the star, the 
solution that remains finite is given by / cx x~^ , which is that of a point dipole. (The case s = 2 corresponds to a constant 
magnetic field in the z direction.) 

Since the density of charged particles decreases to zero at the surface of the star, there cannot be surface currents and 
the current density has to approach zero at the surface, implying that equation (I18p must be satisfied also at a; = 1. In 
addition, the magnetic field must be continuous across the surface, which implies that both / and /' should be continuous. 
Since / oc x~^ outside, it then follows that 

f = at x = l . (19) 

Since /(I) 7^ 0, this equation requires that |/(a;)| decrease locally towards the surface. Moreover, we have |/(0)| = ^ l,/(l)li 
which implies that |/(a;)|, or equivalently |d(a:,&)|, has at least one maximum somewhere within the star. 

In addition to the boundary conditions at the stellar surface (equations 1181 and I19p . the function / must also satisfy 
regularity conditions at the center. Since the force density must remain finite, both the magnetic field and the current density 
must remain finite as well. In particular, for a trial solution of the form f cc x" , we have iJpoi oc x'''^ {2r cosO ~ sO sinO) 
(equation [15}, and 47rJpoi/c oc — (s + l)(s — 2)x"~'''4>sin9 (equation [16}. In order to avoid singularities and multi-valued 
functions at the origin, we must have either s = 2 (corresponding to the zero current case), or s > 3. Consistent with this, we 
seek a solution of the form 

fix) = f2x' + /4X* + fax" . (20) 

We need at least three terms in this polynomial ansatz, in order to be able to satisfy the two homogeneous boundary conditions 
at the surface (equations 1181 and I19p . Considering the solution outside the star, and normalizing /(I) = 1, we then have 




(21) 

Poloidal field lines are lines of constant a, and are illustrated in Fig. [2] for the field configuration discussed here. Note 
that, for a given x, a is largest along the equator. It increases smoothly from at the center, reaches a maximum at 



X 



max 



(14 — \/2T)/15 ^ 0.792, where its value is dmax = ./max = (931 + 21^21)7900 ^ 1.14, and then decreases back down 
to 1 at the surface. The equatorial circle of radius Xmax is known as the neutral line, and the poloidal magnetic field vanishes 
there. On the other hand, the equation aix, 9) = 1 defines the last magnetic surface that closes within the star. Consequently, 
the region where 1 ^ q ^ Qmax is occupied by field lines closing inside the star. 



2.4 Toroidal field 

An axisymmetric magnetic field with poloidal and toroidal components can be written as (equation [10} 

B = B, (vpoiVa X Vfli + rytor/3V<;!>) , (22) 

where T^poi and »7tor are dimensionless constants that determine the relative strengths of the two components of the magnetic 
field. As discussed in i\2.2l /3 must be expressible as a function of a. Moreover, the toroidal field must vanish outside the star, 
since there are no currents to support it there. Since both a and /3 are constant along the poloidal field lines, it follows that 
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Figure 2. Magnetic field lines for a poloidal field given through equations II14I I and I I21I I. a is constant along the field lines. The toroidal 
field is present only in the shaded donut-shaped region within the star. The stellar surface is shown with a dashed line. The field outside 
the star is that of a dipole (which is curl- free, i.e. there are no currents outside the star). 



the toroidal field is non-zero only in a torus-shaped region defined by the poloidal field lines that close inside the star. The 
boundary conditions on the poloidal field remain unchanged, and we can still use the results of the previous section (equations 
ll4l and [2T1l . In this case, the last poloidal field line that is closed within the star is given by a{x, 9) — n/2) = 1. Consider 
a simple relation of the form 

"^1' (23) 

^ for a < 1 . ^ ' 

In order to avoid fast Ohmic dissipation, the current density inside the star must be continuous across the boundary where the 
toroidal field vanishes. This implies that we must have n > 1, so that the current due to the toroidal field decreases smoothly 
to zero at the boundary. In this paper we will consider the case n = 2, but whenever possible we will keep track of the power 
n for completeness. The equilibrium pressure and density perturbations corresponding to this field structure are calculated in 
Mastrano et al. (2011). 

The extent of the toroidal region in the meridional plane can be calculated from the condition that a{x,0) 1. The 
radial extent is largest along the equator, and its limits are given through the roots of f{x) = 1 in the interval ^ a; ^ 1, 

which are x = (27 — \/249) /30 ~ 0.612 and x — 1. The largest angular extent is given by the condition l//max ^ sin^ 6, 
where /max ~ 1.14 as noted in the previous section, which yields 1.21 < S < 1.93 in radians (or, 69.4° < S < 110.6°). The 
region where the toroidal field is present is depicted in Fig. (2] 



2.5 Amplitude of the magnetic field 

The poloidal field is largest along the axis and has a maximum at the origin, while the toroidal field is largest along the 
equator and has a maximum at a; ~ 0.782. In our notation, the largest amplitudes of the two components are 

35 

(-Sp°l)max = —VpolBo = fepolBo and (Btor)max ~ 0.0254 ?7torBo = 6tor-Bo . (24) 

The coefficients 6poi and 6tor defined in this way are dimensionless. We also note that the surface magnetic field (which is 
entirely poloidal) has a maximal amplitude of 2?7poi-Bo at the poles (where it is radial), and a minimal amplitude of r]po\Bo at 
the equator (where it is tangential to the surface). 

Consider the energies stored in the poloidal and toroidal components of the magnetic field, 

£^poi = / ISpoii'dV = ^Boi?2??poi = 2.77 X IQ-^BlRthl,, , 

■' (25) 



- / \BtorfdV ^ 4.12 X IQ-'^B^Rlritor ~ 6.39 x 10"^Boi?*&L 
)7r J 
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The integration for the poloidal part is carried over all of space, while the volume where the toroidal field is present is much 
smaller. 



3 STABILITY 

To study the stability of the magnetic field, consider small fluid displacements around the equilibrium given by equation 

~Q^^ Qoj\ = 5{VP + qV<^ - = -T{£.) . (26) 

Here, 5 denotes Eulerian perturbations due to the displacement field ^, and T is the net force density induced by the 
displacements. Note that there are two types of perturbations in our treatment: the magnetically induced ones with respect 
to the non-magnetic equilibrium, which we denote by the subscript 1 as in equation Q, and those induced by the small 
displacement ^ with respect to the magnetic equilibrium. The latter can be described either as Eulerian perturbations 5 
(changes at fixed locations) or Lagrangian perturbations A (changes as a fiuid element is displaced), which are related 
through A = 5 + ^ ■ V (Friedman & Schutz 1978). 

There are two ways along which one can proceed from equation (|26|) in order to determine the stability of the magnetic 
field configuration. One method is to solve the equation for the perturbations explicitly to determine the frequencies u}, and 
require them to be all real, uj^ ^ 0. Another method is to employ the energy principle of Bernstein et al. (1958), which 
has the advantage that one does not need to actually solve the equation; however, it also has the drawback that it is often 
quite complicated to draw general conclusions. The energy of the perturbations can be written as the sum of hydrostatic and 
magnetic terms, 5W = ~\ j J^dV = SWi-^yd + 5W^mag, where (Akgiin & Wasserman 2008) 

SW^yi =^ J [rP(V ■ if + ii ■ VP)(V ■ ^) - (^ • V$)(V • gi) + • V5$] dV 

-^f[rPV-i + i-VP]i-dS, (27) 



5Wmag =^ 



\5B\^ J-SBx£, 



dV + -^ i [i{B ■ SB) - B{£, ■ SB)] ■ dS . 

OTV 



An c 

The magnetic field perturbation follows from Faraday's law of induction, 

SB = V X i^x B) . (28) 



3.1 Implications of non-barotropy 

The pressure in a non-barotropic fiuid can be written as P{g,s), where s is the specific entropy or chemical composition, 
depending on the type of star (as discussed in 32]). In the non-magnetic background equilibrium, this quantity is a function 
of density, so(f?o), because both so and go are functions of radius. Thus, the background equilibrium is described by a single 
index, 

d In £10 \ din g ) ^ \ dins J ^ dingo 

For the non-magnetic equilibrium described in H2.ll we have 7(2;) = (5 — 3x^)/(2 — x^), which decreases monotonically from 
7(0) = 5/2 to 7(1) = 2. 

For long equilibration times, the quantity s of a given fiuid element remains constant as it is displaced, therefore As = 0. 
Then, the Lagrangian perturbations of pressure and density are related through 

APf^lnP^Ag^^Ag 
P \dlngj^ g g 

Similarly, working to lowest order in (dropping terms of the order ^B^), and using A^i — — £)V ■ i, Sg = —V ■ (g^), 
Ss = — ^ • Vs ~ —{dso/dgQ)i ■ V^po, and the definitions of 7 and F, the Eulerian perturbation of pressure can be written as 

— ~7— + (r-7 — • (31) 

. s g g 

In a non-barotropic fiuid, F 7^ 7, and the hydrostatic force (which we define as the sum of pressure and gravitational forces, 
fhyd ~ ~ i?V<I>) now gives rise to an additional term due to buoyancy, which is proportional to the difference between 

the indices. Upon the application of small perturbations we have 

Sfbyd = - (5(?V$ - g'VS<i> = -g'V + 5$^ + (^^ - 1^ AgV<i> . (32) 

In a stably stratified star F > 7, and the second term acts as a restoring force. Typically, in upper main-sequence stars 
F/7 — 1 ~ 1/4, in white dwarfs F/7 — 1 ~ r7/500, where T7 is the internal temperature in units of lO'^ K, and in neutron 
stars F/7 — 1 ~ few % (Reisenegger 2009 and references therein). 



SP 


/91nP^ 


) ^ + 


/91nP\ 


T ^ 
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3.2 Implications of stable stratification 



Consider the integrands of the hydrostatic and magnetic parts given by equation (|27|) . For simplicity, we will always be 
concerned with cases where the surface integrals vanish (i.e. ^ = at the surface). Moreover, we will employ the Cowling ap- 
proximation of neglecting perturbations of the gravitational potential (5$ = 0), which simplifies the calculations considerably. 
We have 

Shyd = rP(V ■ $f + ■ VP)(V V$)(V • g^) , 

(33) 



^mag — 



— [|5BP -£ X (V X B) ■ 



We can write the adiabatic index of the perturbations as F = Fo 

,'dlnP^ 

i 



Fi, where, from equation ((3^ 



and Fi is the magnetic correction. Note that |Fi|/Fo ~ 
Thus, the hydrostatic integrand can be rewritten as 

£hyd =(Fo - -f)Po{V . £)2 + If (V ■ Qoif 



ding 
Pi\/Po 



(34) 



\gi\/go ~ l$il/"l>o ~ B^/Po < 10"" (Reisenegger 2009). 



(35) 



+ (FiPo + FoPi)(V • if + ■ VPi)(V • - • V$o)(V ■ gii) - (£ ■ V$i)(V • go^) 



For stably stratified stars, Fq > 7, so that the first two terms of the integrand are positive definite. The remaining terms of the 
integrand are corrections due to the magnetic field. These, as well as £mag, can be positive or negative, but their magnitude 
is ^ ^^B^ /L^, where L is some length scale characterizing the spatial variations of the magnetic field. In order for the total 
energy to be negative, thus allowing for the existence of instabilities, the first two terms of £hyd must also be small, 

^2^2 



(Fo-7)Po(V-£)^< 



and 



(36) 



These are constraints that need to be satisfied by the displacement field in order to potentially lead to instabilities. They also 
imply the following bounds for the remaining terms in fhyd, 



|Fi|Po(V-0 -TolPiKV-^)" < 
1(£ ■ VPi)(V • £)| ~ |(£ ■ V$o)(V • eiOl < 



(Fo - 7)^'oi' 



^2^3 



v/(Fo - 7)PoL2 



(37) 



|(£- V$i)(V-po01 < 

Here, we assume that both 7 and Fo are of order unity. Although Fq — 7 ~ 10""^ ^ 1 in some realistic cases, it is still much 
larger than the ratio of magnetic pressure to background pressure, B^/Po ~ 10"". Thus, we conclude that (i) corrections to 
the equilibrium pressure and density due to the magnetic field give rise to terms in the hydrostatic energy that are at least 
a factor of B/ \/ (Fo — 7)Po ^ 10"^ smaller than the (potentially destabilizing) magnetic energy contributions, and therefore 
can be left out; (ii) the conditions given by equation (I36|l also imply that the radial component of the displacement field is 
smaU, Cr/C' < BV(Fo - 7);'o < F 



3.3 Energy of perturbations for a general displacement field 

In this section, we will write down the energy of arbitrary perturbations for poloidal and toroidal fields. Since we assume 
axisymmetry, and none of the equilibrium quantities depends on the azimuthal angle 0, we can express the displacement field 
in general as a superposition of components of the form 

£ = [P(r-, e)r + S{r, 6)0 + iT{r, 6)^ r sin Oe'""'' , (38) 

which can be analyzed separately for different m (as they do not mix in the energy). For notational convenience, we have 
explicitly written out a factor of cylindrical radius. In general, the dimensionless functions R, S and T will be complex, but 
only the real part of £ is physically relevant. Therefore, products should be treated as ZZ*, where * denotes the complex 
conjugate^ 

^ Caution must be taken in using complex notation to describe real physical quantities. Here, we are dealing with functions of the form 
/ = F{r,6)e'""'f' and g = G(r, 6)e*'"'^, and are interested in integrals of the products of their real parts (denoted by R), which can be 
written as 

n{m{g)d<l> = I ma*)d<l, = 7r5R(FG*) . 
Jo ^ Jo 



© 0000 RAS, MNRAS 000, 000-000 



Stability of magnetic fields in non-barotropic stars 9 



The energy of the perturbations can be calculated from equation ((33]). For notational convenience, define an operator A 
and an auxiliary quantity Dm by 

drir-'R) dg{S sin^e) mT 



A{u) = RdrU + 



Sd, 



and 



(39) 



A is the directional derivative along the displacement field, ^ • Vii = A(it)r sinSe""'*, where u is an equilibrium quantity 
independent of the angle 0. Dm is the divergence of the displacement field, V • ^ = Dmr sin ^e™**, and we will explicitly keep 
track of its dependence on m. Defining vj — rsinO, the hydrostatic part of the integrand (equation I33|) becomes 



fhyd = -Tu'^HlrPDrnDm* + [A(P) - gA{$)]Dm'' - A{g) A* {<1>)} 



(40) 



The factor 1 /2 arises as a consequence of the complex notation, as discussed in footnote [T] The equation of hydrostatic 
equilibrium (equation |3J relates the pressure, density, and gravitational potential to the magnetic field. Therefore, Shyd 
depends implicitly on the functions a and /? for the poloidal and toroidal components of the magnetic field through the small 
corrections that these induce on the background quantities. However, as discussed in i\3.2\ these corrections can be dropped in 
the calculation of £hyd, so the equilibrium quantities P, g, and <1> in equation (|4Up can be taken as their non-magnetic versions 
Po, £'0, and $o. 

The calculation of the magnetic part of the integrand is more involved. The perturbations of the poloidal and toroidal 
components of the magnetic field (equation I10|l . are given by equation (|28|l as 



SBpoi =V X (^ X Bpoi) = V(^ ■ V0) X Va - V(^ ■ Vq) x 

( mTdoa — dg[vjA{a)] ^ mTdrCt — dr[TuA{a)] ' i{drTdga ~ deTdrO) -i\ imtf, 



1 



ruj zu r 

SB,,,. = V X (^ X Bto,) = ipVcb • V)^ - (/JV0)(V V)(/?V0) 

imRp^ , imSP^ dr(rRl3) + de{SP)j\ ,m4> 
V H f7 — (p > e 

VO VJ ^ J 

Also, using equation we have 

^ X (V X Bpoi) = ^ X (-AaV(?!>) = -Aa{Sf - RO)e^""'' , 

$ X (V X Bto.) = $ X (V/3 X V0) = ■ V^)VP - ■ VP)V<f> = [iTV/? - A(/?)<^]e™^ 



(41) 



(42) 



The magnetic part of the integrand (equation 1331) can be written as a sum of three terms: one that is entirely due to 
the poloidal field, one entirely due to the toroidal field, and a third term that is a combination of the two components, 

where 



•Spoi = ['5Spoi • <5S;„i - ^ X (V X Bpol) • 5b;oi] 
£tor = [SB,,, ■ SBl,, - ^ X (V X Bto,) • SBl,,] 



(43) 



cross 

After some algebra, we obtain 

^poi 



— K [<5Bpoi • SBl,, + SB,,, ■ SB;,i - ^ X (V X Bp„i) • SB*,,, - ^ x (V x B„, 



sb;,^] 



1 

'8^ 



mTdrO - dr\vjA{a)\ ^ RAa 



+ 



mTdea - 9s[t*7A(a)] _^ SAa 



+ 



drTdga - dgTdrO 



{\R\^ + \S\^){Aaf 



P[dr{rR) + deS] A(/3) 
r 2 

iT 



mPT A(/?) 



+ 



m^PWRC + \S\'' + \T\^) 



(44) 



— 5R I ^ [dr[vjA{a)]del3 - de[vj A{oL)]dr p\ * + ^ [dr{rR) + deS\ \drTdea - deTdrC^ * 
2iml3 



+ 



D c* A , mTdra - dr [vj A{a)] mTdga - dg [wA{a)] 



3.4 Stability of a toroidal field 

In this section, following the derivation of Tayler (1973), we consider the problem of constructing a displacement field that 
makes a purely toroidal magnetic field unstable. In other words, we want to find ^ for which fhyd + £tor < 0. Then, in the 
following section we will examine the stability of the poloidal part for the same displacement. 

As demonstrated by Tayler (1973) for the purely toroidal field, the real and imaginary parts in the energy separate into 
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two equivalent terms. Consequently, it is sufficient to consider the case of real R, S, and T. The function T appears only 
algebraically in the hydrostatic and toroidal parts of the integrand. We then have, from equations (|40|l and (144 |l . 



where we define. 



o 



= 2 [E2{mT)' 



EimT + Eo 



(45) 



E2 =rp , 

Ei=- 2vjTPDo 



iz'A(P) + vjqK{^) 



PHP) 



En 



2 1 



vj TP + — 
47r 



Dl + 



=A(p)A($) - 



m^A{P) - zu^gA{<S>) + 
l3A{/3)A{zu 



I3A{I3) P^Aizj) 



2ttzu 



47r TTVJ 

■Kzo'2 Atvzj^ 



Do 



(46) 



Here, Do = -Dm + mT/rsin9 (equation I39p is the only term that contains derivatives of the functions R and S. The above 
terms can be somewhat simplified using the equation of equilibrium for purely toroidal fields, which follows from equations 
(O, (Hg, and (O as 



A(P) + pA($) = 



PMP) 



(47) 



Since E2 > 0, the integrand given by equation (|45p can be minimized with respect to T for m 7^ 0. In the minimization, 
we hold 7? and S (and therefore Do) constant. The minimizing value is mT/vo = ~Ei/2vjE2 = -Do — gA{^)/TP and the 
minimum of the integrand is fhyd + Star = Eo/2 — -E1/8-E2. Using equation (I39p . this minimization corresponds to setting 
Dm = £'A("l?)/rP, which can be alternatively expressed as 



rpv ■ ^ = f?^ ■ V$ . 

Dropping magnetic corrections to the background quantities (which give rise to terms of the order ^B^), using Sg — 
Aq — — f?V ■ ^, and equation (|3ip . this can be rewritten as 



5P 
P 



. Sg Ag 
Q Q 







(48) 
-V-(eO, 

(49) 



Note that Ehyd is a quadratic function of T (equation I40|) and ftor is a linear function of T (equation [44} . This implies that 
both fhyd and Ehyd + iStor can be minimized with respect to T for m 7^ 0. In fact, the minima of the non-magnetic case (which 
corresponds to minimizing £hyd) and the purely toroidal case (which corresponds to minimizing £hyd +ftor) are both obtained 
for the condition given by equation H48|) . These minima are not precisely identical since the background quantities differ by 
a small amount between the two cases. Equation (|49p implies that the minimum is obtained by setting 5P = 0, which in a 
barotropic fiuid (F = 7) further implies that 5g = 0. The minimum of fhyd to lowest order is fhyd ~ (I/7 — l/r)(^ ■ VPo)^/Po 
(equation ISSp . which is zero for a barotropic fluid, while for a stably stratified non-barotropic fluid it is positive (as long as 
R^Q). 

On the other hand, for m = 0, we have £hyd + i^tor = Eo/2. Thus, in general, we can combine the two cases (m = and 
m 7^ 0) and write the energy for any m as 



^hyd ~r t^tor — ^ 



(l-'^™o)^ 



where 



for 
for 



m = 
m / 



(50) 



We can further rewrite the integrand by grouping the Do terms together and writing them as a complete square, thus 
separating the derivatives of R and S and leaving out only algebraic terms. Defining Km ~ SmoTP + 13^ /Atvzu'^, we have 



--hyd 



1 2^ 

— ZU K-rr 

2 



Do 



1 



5mOQA{$) + 



l3^A{vj) 



2'KVj'^ 



+ ^Uj'^{amR^ + bmRS + CmS'^) 



(51) 



Keep in mind that this integrand is already minimized with respect to T for m 7^ 0. The first term is always positive, and 
the second term forms a quadratic in R and S. The positive definite term can always be made to vanish by a suitable choice 
of the displacement field. Therefore, the integrand is always positive if the quadratic is positive, which corresponds to the 
conditions 



am > 



> 



and 



bm < 4am c„ 



(52) 



These are sufficient and necessary conditions for the stability of the toroidal field (Tayler 1973). Note that they are not 
independent: one of the first two, together with the last one, imply the remaining condition. The coefficients for any m are 
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given through 



drgdr^ - (1 - 5„0) -TT^ 5mO£'9r$ ' 



TP K,„ V 27rr3sin^6» 



2Tvr^sm^e Trr" sin^ 61 47rr4 sin'' 61 ' 

^™ = i^-Smo) — — 

' ^Sr^OOdr^+ ^j' ^ (S^0gds<l>+ Sj"-1. ) (53) 



rKm \ 2nr^ sin 9 J \ sin 

I3drl3 cosO pdeP 2/3^ cos 9 



27rr3 sin^ 6* 27rr4 sin^ ( 
dggde^ , .g'ide^f 1 f, /^^cosS 



+ — r--rT + 



27rr4sin^6l vrr* sin* 61 47rr'* sin'' 61 ' 

These are equivalent to the results given by Tayler (1973), Goossens & Veugelen (1978), and Akgiin & Wasserman (2008), 
albeit the notation is somewhat different. (Here, we have combined the cases m = and m 7^ into a single general form.) 
We have \gdr^\ ^ \^drg\ ^ Po/Ri, and \gdg^\ ~ \^dgg\ ~ B^, so that, to leading order, the coefficients are 

1 1\ {drPof _ „__,,2 (r-7)P0 „ , ,„ , 2 



^ I D - £'oA*' ~ 52 ^^'^ I'*™! ' ~ = t'o'^A • (54) 

Here, A'^ is the Brunt- Vaisala frequency, and lja is the Alfven frequency (i.e. the inverse of the Alfven crossing-time for the 
star). In a stably stratified star, F > 7, where the two gammas, defined by equations (|29p and (|30p . are of order unity. In this 
case, the condition a„i > is comfortably satisfied, and the problem reduces to showing whether 6^ < 4am Cm is satisfied, 
since, if it is true, then the remaining condition Cm > follows trivially. However, note that, when Cm < 0, the magnetic field 
is always unstable, immaterial of the value of bm- The field can also be unstable when Cm is positive, but sufficiently close 
to zero, while bm is sufficiently large (0 < Cm < b'^/iam ~ B'^ /6iiY^ (T — ^)PoR^). This is a very narrow interval. For larger, 

< 4a,„c,„ will always be satisfied. 



3.4-1 Coefficients for m = 

For future reference, we quote the coefficients for m = here. In this case, we haveifo =VP + I3'^/A TTUJ , and the coefficients 
given by equation (|53p reduce to (Tayler 1973; Goossens & Veugelen 1978; Akgiin & Wasserman 2008) 



Ko \ 2nr-^ sm 9 I 2-Kr^ sm 9 



6o = ^^-^-^U$ + ^^lU<^+ /3^cos. 



r 



rKa V 27rr3 sin^ 9 j \ 2Tvr^ sin^ 9 



Pdr 13 cos 9 pdeP _^ 2^^cos^ 



(55) 



sin^e 27rr4sin2 6' 7rr4sin='6l ' 
^„^^g^ 1 ( 13^ cos 9 Y l3dgl3cos9 fi^ cos'^ 

7^ T^A^" 27rr2 sin» 9 ) ^ 2-Kr^ sin^ 9 -kt^ sim 



3.4-2 Coefficients for m ^ 

For m 7^ 0, the integrand given by equation (|5H) reduces to 



/^^ 2 1 2 2 2 

Shyd + Star = TT [dr{rR) + Og S] + -TU {am^oR + bmjtoRS + CmjtoS ) . (56) 

Sivr'' 2 



The coefficients are given through (Tayler 1973; Goossens & Veugelen 1978; Akgiin & Wasserman 2008) 

a,„^0 = -drgdr^ - Vo - .Jl.l2a + 



VP 27rr3 sin^ 9 47rr4 sin" f 

^ _ drgdg^ dggdr^ 2g'^dr^dg^ pdrf3cos9 pdgfi /j-^n 

^ r TPr 27rr3 sin^ 9 ~ 2-Kr^ sin^ a ' 

e^(9fl$)^ l3del3cos9 m^p'^ 

Cm^o - -2 jTp^a 2nr* sin^ 61 47rr4 sin* 
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3.4-3 Proof that all continuous toroidal fields are unstable 

Tayler (1973) shows that an instabihty exists regardless of field strength, if the field has the right structure, and that a toroidal 
field with a non-zero current density on the axis is necessarily unstable. Goossens, Biront & Tayler (1981) further show that 
a toroidal field is unstable if there is some point in the star where the field strength is zero {B — 0), but the derivative of B 
with respect to sin6' is positive. Next, we show more generally that, in fact, all physically relevant toroidal fields are unstable. 

A toroidal field is unstable if, for some value of m, c„i < somewhere in the star. For m = (equation [SSJ , neglecting 
perturbations of the order of B^ caused by the magnetic field, we have 

PdeP cos e p^cos'^e _ sing cos 6*^ f 
sin^ e 7rr-i sin" 9 47rr'* 



On the other hand, for m 7^ 0, we have, from equation (|57[) . 

l3dg/3 cos e _ ta.n"'^-^ edej/S'^ cot""" 6) 

Wo - ~ 27rr4 sin-'* 6 ^ Atvt^ sin^ ^ Atvt^ sin^ 9 



(59) 

Considering specifically m = 1, this simplifies to 

deiP^^ (60) 
47rr4sin^e» ^ ^ 

Now, we need to look at the behavior of /?. In this case, the magnetic field and current density are B = /3V(j!) and 
A-kJ/c = V/? X Vi;/), respectively (from equations llOl and llH) . At the axis, the magnetic field and current density cannot have 
and 4> components, which implies that /? must go to zero faster than zu oc sin 61. Then, as tij — >■ 0, we have A-kJ jc — > ^^^9^/35. 
If we want the latter to be finite, we need /3 to go to zero at least as fast as vj^ oc sin'^ Q as we approach the axis. Thus, in the 
above equations (for m = and m 7^ 0), easily cancels the singularities due to sinS at S = and tt. 

In particular, consider the coefficient co as given by equation (|58|) . If /3 oc sin'^ Q then co vanishes everywhere, which, in 
the best case, implies marginal stability (if &o = as well) . If, on the other hand, /? goes to zero faster than sin^ Q (as is the 
case for the field considered in this paper, for which /? = identically in a finite range of 0), then the function fP" j sin'' 
increases from zero to a finite value somewhere in the interval < 6* < 7r/2 (i.e. it has a positive derivative while cosS > 0), 
and decreases from some finite value to zero for 7r/2 < 6 < tt (i.e. it has a negative derivative while coa9 < 0). Therefore, in 
both cases we will have some regions where co < 0, thus leading to instability. 

On the other hand, as can be seen from equation (|60p . we will have ci < whenever the derivative dg{/3^ cot 6) is positive. 
Note that cot 6^ = at 6I = 0, tt/2, and tt; cot 6^ ^ for < 6* < tt/2, and p'^cotO < for -n /2 < 9 < tt. Thus, if /3 is 
a continuous function, we will have ci < somewhere in whichever hemisphere /3 has some non-zero values. In other words, 
the m = 1 instabilities will happen in those regions where, moving on a spherical shell of constant radius in the direction of 
increasing 9, cot 9 increases from zero to some finite value while cot > 0, or from some finite negative value back to zero 
while cot 6* < 0. 



3.4-4 Application to our particular magnetic field structure 



Now consider the application to our choice of toroidal magnetic field, where /? is given by equation (|23|l . and a is given by 
equations p4|l and (|2ip . Keep in mind the renormalization of the functions a and j3 carried out in accordance with equation 
(|22p (we are now considering the case rj^oi = and ?7tor = 1, so we have only a toroidal field, despite taking a 7^ 0). For 
completeness, in the next few lines, we will keep track of the power n defined in equation (|23|l . Thus, oi{x,9) = f(x)sm^9 
and P{x,9) = [a{x,9) — 1]" in the region where the toroidal field is non-zero. Note that this /3 vanishes long before reaching 
the axis. We consider only the region where /3 7^ 0, i.e. a > 1, or equivalently, l/f{x) < sin^ 9^1. From equations ((58} and 
(|6Up . we have 



co.i 



'■9-1]^ 



Bl/AnRl 



-Go,i 



where Go{x,9) = Acos'^ 9[{n- l)f{x) sm'^9 + 1] and Gi{x,9) = f{x)sin^9{'ir 



(61) 

' 9 ~1) + 1- In equation (|6ip . the coefficient of 
Go,i is always non-positive, and the instability condition co,i < thus requires Go,i > 0. Note that for our particular choice 
of the magnetic field Go ^ 0, thus Co ^5 0, i.e. there is always an m = instability. For m — 1, the instability condition Gi > 
is satisfied for sin^ S « l/f{x), but is not satisfied for sin^ S « 1, i.e. ci changes sign somewhere between these two values. 
The point where the sign change takes place is given as the real root of Gi — 0, 



1 



1 

4n 



1 

4n 



+ 



(62) 



1.14 



The unstable region is l/f{x) < s'm^ 9 < sin^ ^c, and the stable region is sin^ 9c < sin^ 9^1- For n = 2 and f{x) = /m 
f H2.3p . corresponding to the largest extent in colatitude, the stable region is 83.2° < 9 < 96.8°, and the rest of the interval 
where the toroidal field is present, 69.4° < 9 < 110.6° f i]2.4p . is unstable. The contours of the coefficients Co and ci are shown 
in Fig. El 
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Figure 3. Contours of Tayler's coefficients cq (equation I55I I and ci (equation I57I I. The coefficients are sfiown in units of b^^^ / 4it . 
The equator (horizontal) and axis (vertical) are shown in units of stellar radii {vu/Ri, and z/R-t,, respectively). The outer boundary is 
defined by the last poloidal field line that closes within the star (and is tangential to the surface at the equator); both /3 = and = 
on this boundary. This defines the region where the toroidal field exists. The toroidal field is of the form given by equation I I23I1 with 
n = 2. The contours are shown in the range —1.2 (dark) to —0.3 (light) for cq, and —0.9 (darkest) to 1.2 (white) for ci in increments 
of 0.3. Cm < implies instability, cq is zero along the equator (9 = 7r/2) and on the boundary, and is negative everywhere else. The 
maximum of ci occurs along the equator, at a; 0.772; its minima are sA, x ^ 0.772 and 9 ^ 11/2 it 0.205; its radial extent along the 
equator is 0.612 < a: ^ 1; and its angular extent in the meridional plane is 1.21 < S < 1.93. 



3.4-5 Proof that the limiting case of perfect stable stratification implies stability 

Typically, the hydrostatic force is much stronger than the magnetic force. This impUes that, in a stably stratified star, 
any radial displacement will be acting against a prohibitively large buoyancy force. In the limiting case of perfect stable 
stratification, let's consider a displacement field that is perpendicular to the restoring hydrostatic force, which to lowest order 
points in the radial direction (equation [32]). In addition, we require the fiuid to be incompressible. In other words, the density 
remains constant as a fluid element is displaced (i.e. = — £)V ■ ^ = 0, while AP 7^ 0, implying that T — > 00 from equation 
l3Up . These two conditions, namely incompressibility (V • ^ = 0) and orthogonality to the radial direction (r ■ = Q), imply 
that the displacement field is described by a single unknown function (instead of three, as is the case for an unrestricted 
vector field). From these assumptions it also follows that the hydrostatic part of the energy vanishes (to first order in , 
as in equation I35p . Therefore, all we are left with is the variation in the magnetic energy. In our notation, the requirements 
V • ^ = and r • ^ = correspond to setting Dm = and _R = in the integrand of equation (|45[) . To first order in , the 
integrand reduces to (for any m) 



^hyd ~\~ ^tor 



Snr^ 



(63) 



This displacement field is restricted to such an extent that it is not possible to make the first positive definite term vanish by 
a suitable choice, unlike in the general case. Therefore, the coefficient of is no longer sufficient to assess stability. In fact, 
it is possible to show that the integral is always positive, which is not immediately obvious from the above form. It is made 
clearer by rewriting the integrand as 



^hyd "t £tor — 



Sivr^ 



(m^ - 1)3^/3'^ 



+ fidgS + Scote) 



2 deiS^P'^ cose) 



(64) 



The last term integrates to zero, since 5* sin 6* — >■ and f}/ sinO — >■ on the symmetry axis. On the other hand, the sum of the 
first two terms is always positive for ^ 1. For m = 0, the only displacement field consistent with V ■ ^ = and r • ^ = 
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that does not diverge is of the form ^ = £,<ii(i>, which has no effect on the toroidal field (equation [28]). Thus, we conclude that 
the equilibrium is always stable to this very restricted set of perturbations, as previously noted by Dicke (1979). Therefore, 
in order to obtain instabilities, the restrictions due to perfect stable stratification must be relaxed. 



3.4.6 Constraints on the destabilizing displacement field 

In this section we will consider the properties of the displacement field that destabilizes the toroidal magnetic field. Some 
simple observations can be inferred by noting that the energy integrand can be written as the sum of a positive definite term 
and a quadratic as in equation (jSip . We would like the quadratic to be negative and the positive definite term to be as small 
as possible. 

First, consider the quadratic, Q{R,S) = UmR^ + bmRS + CmS^ . Since am is always large and positive, and Cm is small 
and negative in some region, \R\ must be small compared to l^l in order to allow the energy to become negative. On the other 
hand, R cannot be zero, as that reduces the integrand in equation (|56p to the form given by equation (|63p . which was shown 
to be always positive. This is because the integrand in equation (|56|) is already minimized with respect to T, i.e. we have 
implicitly substituted the condition given by equation (|48p . which implies that if i? = then Dm ~ 0, thus leading us back 
to equation (|63|) for the perfect stable stratification. Thus, |-R| must be small compared to \S\, but non-zero. Since Om > 0, 
Q{R, S) can be minimized with respect to R. The minimum is given by 

-Rmiii = -TT^-S* and Qmin = f --r^ + Cm]S'^ . (65) 
2am V 4am / 

For instability, we must have Qmin < 0, implying that bm > 4,amCm. Since typically \bm\ ^ \cm\ dm (equation I54|) . this 
will be satisfied when < (plus a thin region where Cm but very small), and we will have Qmin ~ CmS^ . \i Cm < Q is 
confined to a region of size (Ar, A6') in the relevant coordinates, the displacement field should also be roughly confined to this 
region, as there would otherwise be a positive contribution to the energy from the region where Cm > 0. Thus, in particular, 
\deS\ > l^l/AS (since S must vanish near the boundary of the region, we have \AS\ ~ \S\). 

The positive definite terms for m = and m 7^ are significantly different, and the two cases need to be treated 
separately. First, we will consider the case m 7^ 0. Since the region we are interested in is near the equator (i.e. sinO ~ 1 
and cos6' ~ AS), we have CmT^o ~ — /?^/47rr* (equation l6Up . For instability, we need Shyd + £toi < in equation (|56p . which 
implies [dr{rR) + dgSf < 5^ Thus, we need \dr{rR) + deS\ < \S\ < \deS\Ae, using the above inequality. In other words, the 
more confined the displacement field is in latitude (i.e. the smaller the range A^, forced by the condition Cmjto < 0), the more 
precisely the two derivatives on the left-hand side need to cancel each other. In the limit A^ — )■ 0, it is necessary to enforce 

drirR) + deS ^ . (66) 

In addition, from the confinement to an interval Ar we have \dr{rR)\ > r\R\/Ar, and from the cancelation with deS we 
have \dr{rR)\ > \S\/Ae. The ratio of these two lower bounds is, using equations JSH and (gSJ, $H = (| J?j/|Sj)(rA6'/Ar) ^ 
[\bm^o\/'2am!/^o){rA9/Ar) ~ {B'^/8TvPo){rA9/Ar) ^ 1. Thus, the lower bound on \dr(rR)\ from the requirement of canceling 
deS (even if not exactly) is much larger than the bound from being confined to the region {Ar,A6). Thus, the length scale 
of variation of R must be Sr <^ Ar. Note that this does not mean that ^ is confined to a region as thin as Sr. It could extend 
over the whole Ar, but it would have to oscillate on a radial length scale Sr. 

For m = the coefficient of the positive definite term in equation (|51|l is much larger than in the m 7^ case, Ko 3> Kmj^^o- 
This implies that the quantity in parentheses must cancel out even more precisely. Keeping only leading order terms, we get 

^ + ^^i^ + ^i^^O. (67) 
r3 r sm'^ 9 TPo 

This is different from equation (|66p for the m 7^ case. In what follows we will consider only the simpler case of m 7^ 0, and 
the m — case will be left for future work. 



3.4.7 Particular displacement field for m 7^ 

As discussed in the previous section, the quadratic part of the integrand for m 7^ (equation I56p can be made negative by 
a suitable choice of the amplitudes of the functions R and S. In addition, the first term, which is positive definite, can be 
minimized by a suitable choice of the derivatives of these functions. In particular, the best choice might be when this term 
is made to vanish (Goossens & Tayler 1980; Goossens & Biront 1980), which leads to the condition given by equation H66p . 
This equation is satisfied by solutions of the form 

R = ^ and S = -djl , (68) 

X 

where x is the dimensionless radial coordinate x — r/Ri,, and \l{x,9) is some scalar generating function for the displacement 
field. While we can choose R and S so that the positive definite term vanishes, T has a particular value for which the integrand 
in equation (|45p is minimized. This value of T is expressible in terms of R, S, and their derivatives, corresponding to the 
condition given by equation (|48p . Using equation equations p9p and (|66p . and keeping only the lowest order terms, we have 

Z^-^o-4P^f^l^l^+f^V. (69) 



TP \r FPo 
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Figure 4. Contours of a greatly exaggerated sample generating function 11 for the displacement field (equation [72]l. The contours are the 
streamlines of the displacement field. The stellar surface is shown as a solid line. The values of the various parameters used in this plot 
are xq = 3/4, Bq = n/3, 5r = 1/5, Sg = 7r/5, and it = 3. The actual generating function used in our calculations is much more confined. 



As noted in i]3.2l to lowest order, we can drop all magnetic corrections to the background quantities in the hydrostatic 
part of the energy (equation I40|) . which then becomes, using the equation of equilibrium (equation [T]), the definition of 7 
(equation [29|l. and the value of Dm from equation (|48l) . 

c ^ l^ (c'r-Po)^ r)2 2.2fl 1 r^2 2 ■ 2 a /-rnN 

thyd ~ n -R r sm e ^ -am^oR r s\a 6 . (70) 

2 \7 1 / ft) Z 

The last equality follows from equation (|54|) . In a stably stratified star F > 7, so that the hydrostatic part is always positive. 
The toroidal part of the energy follows from equation (|44p . to leading order and using \R\ <g; IS*], 

£tor ~ ^Cm^oS'^r'^ sin 9 . (71) 

Thus, the total integrand given by equation (|56|l reduces to Shyd + Star ~ ^{a-m^oR^ + Cm#0'S'^)r^ sin^ 6. Since |_R| ^ \S\ and 
Urajto 3> l&mT^oj ~ |Cm/o| (equatiou I54p . it follows that the hm^oRS term in the quadratic can be dropped. 

We can now proceed to construct a particular displacement field that will make the purely toroidal magnetic field unstable. 
We will assume that the displacement field is confined to a region within the star and is zero everywhere else. In order to 
prove that the stability conditions are both sufficient and necessary, Tayler (1973) assumes a particular solution of the form 
n(a;, 0) oc sin kx sin 19 in a finite volume, bounded by a surface on which Yl{x,9) = 0. This corresponds to a finite displacement 
field which is tangential to the boundaries. While this form is acceptable for a purely toroidal field (since both ^ and B are 
tangential to the surface), in our case we will eventually incorporate a poloidal field as well, and any discontinuity in the 
displacement field at the boundaries would cause divergences (cutting the field lines). In order to avoid such pathologies, we 
would therefore like ^ to go to zero at the boundary, and its derivatives to remain finite everywhere. Thus, we choose 

^(x,e) = -|[l-X^(:c,e)]^ where x\x,e) = ^-^^f^ + ^^^^ ■ (72) 

The factor ^o/R* sets the amplitude of the displacement field (equation I38p . The displacement field is zero on the boundary 
(defined by x = 1) and outside of it. This particular choice of x corresponds to a donut-shaped region with a meridional 
cross-section in the shape of a distorted ellipse (Fig. [4]). In order for the derivatives of the displacement field to remain finite, 
we must have cr ^ 2. 

As shown in ii3.4.3l the second condition in equation (|52[) is always violated somewhere for any non-singular toroidal 
magnetic field. Therefore, we will have Cmjto < in some region (Fig. (H)). If we choose the displacement field to be confined 
near the minimum of Cm#o and make Cm^oS'^ to be the dominant term in the quadratic in equation (|56|) . then the energy 
will be negative. Thus, we want, roughly, amjtoR^ < \cm^o\S^ ■ For a displacement field given by equations (|68|) and (|72p . this 
implies that we must have |7?|/|S| ~ Sr/Sg < \/\cm^o\/o.m^o ~ B/y'To ^ 1. Here, we have used \x — xo\ ^ 5r, \9 — 9o\ ^ 5g, 
and equation (|54p . Since 5g cannot be much larger than the angular extent of the negative region of the coefficient Cm/o (Fig. 
[Sjl, this then imposes a very stringent upper limit on Sr. 
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3.5 Stability of a toroidal field in the presence of a weaker poloidal component 

The displacement field constructed in the previous section makes the sum of the hydrostatic and toroidal parts of the energy 
negative, thus leading to an instability. In the present section we will consider the case when a weaker poloidal component is 
added. This poloidal field will give an additional positive contribution to the energy and will help stabilize the instability of 
the toroidal field. Our goal is to determine the minimum strength of the poloidal field relative to the toroidal field in order 
to achieve stability. Note that this treatment inherently relies on the implicit assumption (validated by the eventual results) 
that the poloidal field is sufficiently weaker than the toroidal field, so that the displacement field discussed in the previous 
section is still close to being the most unstable mode. 



3. 5. 1 Leading order estimates of the energy terms 

Here, we will give estimates of the hydrostatic, toroidal, and poloidal parts of the energy in terms of the parameters of the 
particular displacement field constructed in the previous section. Since our R, S, and T axe real functions, the cross term in 
equation (|44|l has no real part and is physically irrelevant. The total energy in terms of the integrand £ is (from equation (27)), 
carrying out the integration over 0, 

5W = SWhyi + SWtor + SWpol ^^J £dV = n J £r^ sin OdrdO . (73) 

The hydrostatic and toroidal integrands are given to leading order in /8nPo by equations (|70|l and (I7ip . On the other 
hand, the poloidal integrand is given by equation (|44p . Here, we need to make use of the azimuthal displacement, which 
is related to the other two components by equation (|69|l . Since ^ 15*1 (or Sr <^ 5g < 1), to leading order we have 
mT ~ 2S cosO. Thus, R oc 1/Sg, S (X 1/Sr, and T oc 1/5,. • Each subsequent derivative dx of the displacement field brings in 
an additional factor of 1/Sr, and similarly, dg brings in a factor of l/Sg. It then follows that the four terms in the poloidal 
integrand as given by equation (|44|l scale as 1/5^1 ^/SrSg, l/S^, and 1/5^, respectively. We thus conclude that the first and 
third terms in the poloidal integrand are the largest, followed by the second term, while the fourth term, which is also the 
only negative term in the expression, is the smallest. Thus, it becomes obvious that, for the displacement field of the form 
constructed here, the poloidal contribution is overwhelmingly positive. To leading order, keeping only the first and third terms 
of the poloidal integrand (equation 144 1) . we have 

£poi^^^^(l+'-^) . (74) 



Sirr^ 

We define a new coefficient for the poloidal field, in analogy to the coefficients Omj^o s.nd Cm^o, 

= f f^)7l + i^) . (75) 
47r \r3sm6iy \ J 

Thus, the hydrostatic, toroidal, and poloidal energies can be written as (using equations 1701 and 1711 and the above definitions) 

5Whyd j f^mi^oR^r^ sin^ OdrdO , 

SWtor c^^oS^" Sin' OdrdO , (76) 

SWpol J dm^o r^siri' OdrdO . 

These approximations are remarkably accurate for the typical values of the parameters discussed in ii3.4.7l (the errors with 
respect to the exact integrals are of the order of 10~^ or less). 

In the limiting case of a vanishing area of integration (i.e. as 5^ — > and Sg ^ simultaneously), all slowly varying 
quantities can be taken as constant, and can be evaluated at the center of the displacement field (xo, Oo). In other words, the 
integrations can be carried out over the rapidly varying functions of the variables {x — xo)/Sr and (0 — Oo)/Sg, and all other 
slowly varying functions of x and can be taken out of the integrations. We thus have 

7r_/?^ 2, 3 f 2 

SWhyd —^am^o{xo,Oo)xosm Oo / (dgll) dxdO , 



SWtor —^Cmi.Q(xQ.OQ)xQsm Oq J (OxU.) dxdO , (77) 

SWpol ^ —^dm^o{xo,0o)xlsm'^0o / — ^ dxdO . 

The remaining integrations can be carried out in polar coordinates through the substitutions (x — xo)/Sr = xcosi/; and 
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Table 1. Numerical values of the dimensionless coefficients fcjjyd, fctor, and fcpoi, defined through equation I I79I I. for xq = 0.772 and 
60 = 1.37, corresponding to the minimum of the coefficient ci (Fig. [Sj. The background quantities are taken as in i|2.2l with p = 
8itPc/B^ = 10^, and the magnetic field structure is that of ^2.31 and ^2.41 We have also set m = 1 and r/7 — 1 = 1/4. 



Coefficient Value 



fchyd 2.70 

fctor 1.06 

fcpol 0.0832 



{6 — 60) /5r = xsini/'; yielding 



Se J 9 6(a-l/2)(a-3/2) ' 

As noted following equation (|72p . we are interested in the case a ^ 2. With these definitions, the limiting forms of the energy 
perturbations can be written as 

5M^hyd ^ (r/7 - l)P.Coi?*fchyd/.^ , 

SWtor ^ ^" bLfcto./.^ , (79) 

OTT Or 

XU/ V ^o^oR* 1,2 , J 
oWpol -> — OpolKpol^CT^ , 

where fchyd, fctor and fcpoi are numerical constants which are independent of 5r, Se, and a, and whose values are given in Table 
[1] The factor r/7 — 1 (also evaluated at the same point as the coefficients) has also been explicitly written in the hydrostatic 
part in order to keep track of the dependence on stable stratification. The amphtudes of the toroidal and poloidal fields are set 
by &tor and fcpoi as defined through equation (|24|) . In addition, we have explicitly factored out all dimensional quantities ^o, Pc, 
Bo, and i?*. The convergence of the hydrostatic, toroidal and poloidal energies to the limiting values of the approximations 
is demonstrated in Fig. [S] 



3.5.2 Stability cntena 

For stability, we must have 5W = SWhyd + SWtoi + SWpoi > 0. Using equation (|79[l this can be written as, noting that fctor is 
defined as a positive number and dropping common factors, 

(r/7 - l)pki^ydla-^ - btarktorI<Tj- + fcpolfcpol J'cr ^ > . (80) 

Here, we have defined p = 8tvPc/B^. We can rewrite this inequality as a lower bound on the amplitude of the poloidal field 
relative to the toroidal field, 



7 — > (>e-r- 

Otor / Jct 



ktor 









fchyd (r/7 ^ l)p fSr^ 
^pol ^tor 



(81) 



We are interested in finding the minimum poloidal field strength that can stabilize the magnetic field against all possible per- 
turbations. Therefore, we would like to maximize the expression on the right-hand side with respect to the various parameters 
involved {a, 6r and Se). First, the ratio I^/Ja can be maximized with respect to a, which for a ^ 2 gives a = (3 + \/3)/2 ~ 2.37 
and Icr/Ja = 3(2 — \/3) ~ 0.804 (see equation 1 78 |l . Next, the expression in square brackets can be maximized with respect to 
the ratio Sr/5e (or, equivalently, with respect to Sr while keeping Se constant), yielding 

A \ ^ fc 

\ '^tor '^tor (82) 



Se) 2fc^^^ (r/7 - l)p ■ 

Note that for this value we indeed have 5Whyd + SWtoi < 0, as is required for the instability of the purely toroidal field in 
the first place. Setting 6tor = 1 (which corresponds to measuring the poloidal field strength in terms of the toroidal one), 
r/7 - 1 = 1/4, p= 10^ and using the tabulated values of the coefficients from Table [TJ we obtain Sr/Se « 10 . 

The remaining dependence of (6poi/6tor)^ on 5e is monotonically increasing, therefore we need to evaluate the largest 
physically reasonable value of this parameter. For a displacement field centered at the minimum of the coefficient ci, this 
value cannot be much larger than the angular extent of the negative region of the coefficient (Fig. [3| . Otherwise, the toroidal 
energy SWtor could no longer be made negative. In fact, for any sufficiently small 5,. {5,. < 10^^) this condition (namely, 
SWtor < 0) translates into Se < 0.4. In fact, we will take the largest value of this parameter to be about Se ~ 0.24, which also 
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Figure 5. Errors between the limiting values (equation 17911 and the exact values of the hydrostatic, toroidal and poloidal energies, as 
functions of Sq, where error = (limiting value - exact value)/exact value. All parameters, background quantities, and coefficients involved 
in the calculations are taken as in Table [T] In addition, 5^ and a are taken to have the values that minimize the total energy (or, 
equivalently, that maximize the ratio fepoi/btor, as discussed in the text following equation [SlJ, and Sg is allowed to vary. The exact values 
converge to the limiting values as Sg decreases. 



corresponds to where our approximations start to fail (at this point the relative error between the exact and limiting values 
for the toroidal energy reaches 100%, corresponding to a factor of 2 error, Fig. [SJ. 

Plugging equation (|82p into equation (|8ip . and using the values of the parameters and coefficients discussed here and in 
Table [T] while explicitly keeping track of the dependence on btor, r/7 — 1 and p, we obtain 

^pol ^ r2 fetor ^tor ~ t; S v 1 ^tor foo\ 

Since p/&tor 3> 1, this is a very small lower bound on the amplitude of the poloidal field needed to stabilize the toroidal field, 
thus effectively justifying our treatment of the poloidal field as small in comparison to the toroidal field. We can rewrite this 
result in terms of the energies stored in the magnetic and gravitational fields. From equations @ and (|25p . we have 

^poi ^ 4 3 y ^ ^ g 7 X 10-' ^ . (84) 



Etoi \ btoi J Ef^ 

Replacing these in equation (|83|l . we get 

^'poi 3.7 -Btor (85) 

-Btor ~ — 1 -Egrav 

Observations provide us with an upper limit on the poloidal magnetic field strength. The above equation then gives us 
an upper limit on the toroidal field strength. Taking -Bpoi/Sgrav < 10"^, we get -Btor/Sgrav < 5.2 x 10"'^ ^T/j — 1. Thus, 
the maximum toroidal field strength depends on how stably stratified the star is through the factor r/7 — 1. The more 
stably stratified the star, the stronger the maximum toroidal field for a given poloidal field strength. For main-sequence stars, 
r/7— 1 « 1/4 and we obtain -Btor/Sgrav < 2.6x 10"'', while for neutron stars, r/7— 1 ~ 10"' and we have -Etor/i?grav < 5x 10"^. 
This also implies that a significant portion of the magnetic energy may be hidden in the toroidal field, while only the poloidal 
field is observed. 

In particular, we can apply this result to the case of magnetars. For a I.4M0 star with 10 km radius, the gravitational 
energy is -Bgrav ~ 4 x 10^'' erg (equation [9)l . The energy of the poloidal field is i?poi « 2 x 10'** Bis erg (equation [25} , where 
_Bi5 is the magnetic field strength at the equator in units of 10^^ G. (Note that, given the existence of closed poloidal field 
lines, this number is an order of magnitude higher than the most naive estimate obtained by multiplying the energy density 
corresponding to the surface field, /Siv, by the volume of the star.) Using r/7 — 1 ~ 10"', we then obtain an upper limit 
for the energy of the toroidal field as -Etor ^ 5 x 10*^ Bis erg. Note the linear dependence of the maximal toroidal energy 
on the surface magnetic field strength. Soft gamma repeaters (SGRs) release as much as a few lO*** erg energy in a single 
outburst (Mereghetti 2008). If the outbursts are repeated every century or so over a period of 10 millennia, then the total 
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energy release is of the order of a few 10 erg. These could only be explained in terms of magnetars with poloidal fields in 
the excess of 10^^ G. However, inclusion of a toroidal field increases the potential energy available for outbursts, and a lower 
surface magnetic field of the order of lO'^* G would be sufficient to explain the observations. 

A particularly interesting case is that of SGR 0418+5729, where the inferred surface magnetic field strength is just 
below 10^^ G (Rea et al. 2010). Its observed X -ray luminosity is ~ 6.2 x 10^^ erg/s, and its characteristic (spin-down) age 
is ~ 2.4 X 10^ yr. If the object is assumed to maintain the same level of activity throughout its life, then the total energy 
required would be ~ 5 x lO**^ erg. Using the formulae of the preceding paragraph, the energy content of a 10^'^ G poloidal field 
is Epoi ~ 2 X 10** erg, i.e. more than two orders of magnitude less than what is required. On the other hand, the maximum 
allowed toroidal field energy in this case is Etor ^ 5 x 10*^ erg, which would be quite sufficient. 

Our results are in general agreement with the simulations of Braithwaite (2009) for a stably stratified, non-degenerate, 
polytropic fluid star with 7 = 4/3 (i.e. a polytrope of index n = 3, which is a reasonable approximation for an upper 
main-sequence star) and F = 5/3. Braithwaite finds a stability condition of the form aEmm^/U < -Epoi/Smag ^ 0.8, where, 
£mag = Epoi + Etoi is the total magnetic energy and U — -Egrav/2 is the thermal energy (by the virial theorem). His simulations 
yield a ~ 10 for main-sequence stars, where r/7 — 1 ~ 1/4. Since in realistic stars -Emag/f/ <^ 1, this implies that, while the 
poloidal component cannot be substantially stronger than the toroidal, the toroidal component can be much stronger than 
the poloidal. 

For the lower bound on the poloidal field strength, Epoi is small and we have i?mag ~ -Etor. The condition given by 
equation (|85|l can then be written analogously to Braithwaite as Epoi/ Etor ^ 2ai5tor/£'grav, where a ~ 1.8/(F/7 — 1). For 
main-sequence stars (r/7 — 1 ~ 1/4), we obtain a ~ 7.4, which compares well with Braithwaite's result of a ~ 10. On the 
other hand, for neutron stars (where, we take r/7 — 1 ~ 10^^, which is different than the value 1/400 used by Braithwaite), 
we obtain a ~ 200. 

Notwithstanding the remarkable agreement between the two approaches, they are also notably different. First of all, 
while our analytic calculations can deal with arbitrary (but small) ratios of magnetic to gravitational energy, Braithwaite is 
forced to use a specific, and not extremely small value for his simulations {Emug/U — 1/400, compared to < 10~® in real 
stars). This means that the magnetic force is significantly stronger, and as a result, stable stratification plays a smaller role 
in stability. This also makes the length scale ratios in the displacement field much less extreme (see equation I82p. and allows 
the unstable wavelengths to be resolved, even if only barely. Thus, effectively, the same instabilities should be manifested in 
both treatments. The proportionality of his final result to general values of the ratio -Bmag/C/ is then stipulated. Secondly, 
Braithwaite's grid of values for -Bpoi/-Emag has only four values for each configuration, which differ from each other by almost 
a factor of 2. Therefore, the value of the coefficient a is not much more precise than that. Thirdly, Braithwaite explicitly 
states that the coefficient Cm (for both m — and m = 1) is always positive for the magnetic fields he considers, and that 
the instability of the toroidal field results entirely due to the failure of the condition 6^ < AamCm f equation I52|l . The first 
statement seems to expressly contradict our conclusion that for any realistic toroidal field, Cm becomes negative at least in 
some regions, and is the leading source of instability, as pointed out previously by Goossens (1980) and Goossens et al. (1981) 
for various special field configurations. Fourthly, we do not consider the m = case in this paper. Tayler (1973) conjectures 
that the m = 1 perturbations "seem likely to be the worst instabilities in the linear regime", and Spruit (1999) notes that 
the m = 1 mode "occurs under the widest range of conditions" . However, as Braithwaite reports, there are instabilities that 
arise from the m = mode as well, and need to be considered. Thus, in our treatment we have only one mode, whereas the 
numerical simulations in principle have all modes, so the simulations should be more unstable. It is also possible that, with 
the addition of the poloidal field, the particular displacement field constructed in this paper no longer corresponds to the 
most unstable mode, which would also imply that we are overestimating the stability. Finally, the hydrostatic background 
and magnetic field structure are also not identical between the two cases. 



4 CONCLUDING REMARKS 

Here, we summarize the main conclusions and discuss further implications of our work. 

Stable stratification has an important influence on stellar magnetic equilibria and their stability, by (a) allowing a much 
larger assortment of possible equilibria, and (b) strongly constraining the displacement fields that might destabilize these 
equilibria. We can easily construct simple analytic models for axially symmetric magnetic fields compatible with hydromagnetic 
equilibria in stably stratified stars, with both poloidal and toroidal components of adjustable strengths, as well as the associated 
pressure, density and gravitational potential perturbations. This makes it possible to directly study their stability. For a weak 
magnetic field (in the sense that the Alfven frequency is much smaller than the Brunt- Vaisala buoyancy frequency) , the terms 
in the energy functional involving fluid perturbations due to the magnetic field are small and can be ignored, which makes 
the algebra much simpler in the cases of poloidal or mixed fields. (For a toroidal field, it does not simplify the algebra much, 
but it simplifies the physical interpretation.) 

There is an important difference between the leading order instabilities of toroidal and poloidal fields. In toroidal fields, 
instabilities result from the slipping of magnetic loops around the magnetic axis (Tayler 1973) . These instabilities are strongly 
restricted by stable stratification to surfaces of constant radius, but are not completely eliminated. As a result, a relatively 
weak poloidal field can be sufficient to stabilize the toroidal field (Spruit 1999). On the other hand, in poloidal fields, the 
instabilities result from the slipping of magnetic loops around the neutral line (Markey & Tayler 1973; Wright 1973). In this 
case, stable stratification is of less help in eliminating instabilities, because, while it restricts radial displacements, it does 
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not help with perturbations that are perpendicular to the radial direction, which are just as easily achievable in this case 
(ignoring curvature effects due to the fact that the neutral line is a circle). Therefore, one might expect that a relatively 
stronger toroidal field would be needed in order to stabilize a poloidal field. These are consistent with the upper and lower 
bounds found by Braithwaite (2009). 

Previous literature (Tayler 1973; Goossens et al. 1981) has given proofs of instability for toroidal fields satisfying special 
criteria. In particular, it is often repeated that toroidal fields are unstable near the axis (Tayler 1973; Spruit 1999). Here, we 
prove that in fact all toroidal fields of any realistic structure, in general (barotropic and non-barotropic) fiuids, are unstable. 
This is true even when the toroidal field is contained in a region far away from the axis, as is the case considered in this 
paper. The instability always happens near the high-latitude limits of the region containing the toroidal field (i.e. farthest 
from the equator and closest to the poles), immaterial of the exact shape of this region. This then allows us to construct a 
displacement field that should be a reasonable approximation for the most unstable mode, compliant with the constraints of 
stable stratification. 

We find that the toroidal field instability considered in this paper is stabilized by a poloidal field that satisfies Epo\/Etoi ^ 
2aEtoi / Egiav (eauation l85|) . where a ~ 1.8/(r/7 — 1). For main-sequence stars, we find that a ~ 7.4, which compares well with 
the factor of a ~ 10 obtained by Braithwaite (2009) through numerical simulations. For neutron stars, we obtain a ~ 200. 
Since observations provide us with an upper limit on the surface poloidal field, this result can then be used to place an upper 
limit on the internal toroidal field. We find that the energy stored in the toroidal field within the star can be significantly larger 
than the total energy of the poloidal field, particularly if the latter is weak. Such strong magnetic fields hidden within stars 
can provide a substantial additional energy budget to power magnetar activity, as well as cause significant stellar distortions 
with implications for precession and emission of gravitational waves (Mastrano et al. 2011). 
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